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The theoretical investigation of spectral functions and pseudogap in 
systems with strongly correlated electrons is discussed, with the emphasis 
on the single-band t-J model as relevant for superconducting cuprates. The 
evidence for the pseudogap features from numerical studies of the model is 
presented. One of the promising methods to study spectral functions is the 
method of equations of motion. The latter can deal systematically with 
the local constraints and projected fermion operators inherent for strongly 
correlated electrons. In the evaluation of the self energy the decoupling 
of spin and single-particle fluctuations is performed. In an undoped anti- 
fcrromagnet the method reproduces the selfconsistent Born approximation 
(SCBA). For finite doping the approximation evolves into a paramagnon 
contribution which retains large incoherent contribution in the hole part. 
On the other hand, the contribution of longer-range spin fluctuations is 
essential for the emergence of the pseudogap. The latter shows up at low 
doping in the effective truncation of the large Fermi surface, reduced elec- 
tron density of states and at the same time reduced quasiparticle density 
of states at the Fermi level. 

PACS numbers: 71.27.+a, 72.15.-v, 71.10.Fd 

1. Introduction 

One of the central questions in the theory of strongly correlated electrons 
is the nature of the ground state and of low energy excitations. Experiments 
in many novel materials with correlated electrons [1] reveal even in the 'nor- 
mal' metallic state striking deviations from the usual Fermi-liquid univer- 
sality as given by the phenomenological Landau theory involving quasipar- 
ticles (QP) as the well defined fermionic excitations even in the presence 
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of Coulomb interactions. The focus has been and still remains on super- 
conducting cuprates where there exists now an abundant and consistent 
experimental evidence for very anomalous low-energy properties [1], besides 
the most evident open question of the origin of the high-T c superconductiv- 
ity. 

In particular, the attention in the last decade has been increasingly de- 
voted to the underdoped cuprates, where experiments reveal characteristic 
'pseudogap' temperatures, which show up crossovers where particular prop- 
erties change quantitatively. As schematically presented in Fig. 1, there 
seems to be an indication for two crossover scales T* [2] and T sg [1]. The 
existence of both is still widely debated, in particular whether both could be 
a manifestation of the same underlying mechanism. Nevertheless we refer 
to them (as usually in the literature) as the (larger) pseudogap scale T* and 
the spin-gap scale T sg for the lower one. 
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Fig. 1. Schematic electronic phase diagram of cuprates 

The T* scale [2] shows up clearly as the maximum of the spin suscep- 
tibility Xo(T = T*) [3]. The in-plane resistivity p(T) is linear p oc T for 
T > T* and decreases more steeply for T < T*. The same T* appears 
in the anomalous T-dependence of the Hall constant Rh(T) for T < T* 
[1]. For the theoretical considerations the most straightforward signature 
of pseudogap is the reduction of the specific heat coefficient 7 = Cy/T be- 
low T ~ T* in the underdoped cuprates [4]. Namely, within normal Fermi 
liquids 7 directly measures the quasiparticle (QP) density of states. The 
spin-gap crossover T sg has been identified in connection with the decrease 
of the NMR relaxation rate 1/Ti for T < T sg [1], related to the reduction 
of low-energy spin excitations. Even more striking is the observation of 
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the leading-edge shift [5] in the angle-resolved photoemission spectroscopy 
(ARPES) measurements at T > T c , a feature interpreted as a d-wave SC 
gap persisting within the normal phase. 

It seems to some extent plausible that the T* crossover is related to the 
onset of short-range antiferromagnetic ( AFM) correlations for T < T* , since 
x(T) in an undoped AFM has a maximum at T* ~ 2 J/3 the temperature 
corresponding to a gradual transition from a disordered paramagnet to the 
one with short-range AFM correlations. In this contribution we concentrate 
on our theoretical results which confirm and explain the existence of the 
pseudogap in the model relevant for cuprates, i.e. the t-J model on a 2D 
square lattice. It should be however pointed out that there there are also 
several alternative theoretical proposals, which do not directly invoke AFM 
spin correlations. 

The single-particle spectral function A(k, u) and its properties are of 
crucial importance, since their full knowledge would essentially clarify most 
open questions concerning the anomalous character of strongly correlated 
electron systems. In recent years there has been an impressive progress in 
ARPES experiments [5, 6, 1] (in particular for cuprate materials) which 
in principle yield a direct information on A(\a, u>). In most investigated 
Bi 2 Sr 2 CaCu202+5 (BSCCO) [5] ARPES shows quite a well defined large 
Fermi surface (FS) in the overdoped and optimally doped samples at T > T c , 
whereby the low-energy behavior with increasing doping in the overdoped 
regime qualitatively approaches (but does not in fact reach) that of the nor- 
mal Fermi-liquid with underdamped quasiparticle (QP) excitations. On the 
other hand, in the underdoped BSCCO QP dispersing through the Fermi 
surface (FS) are resolved by ARPES only in parts of the large FS, in par- 
ticular along the nodal (0, 0)-(7r, it) direction [7, 5], indicating that the rest 
of the large FS is truncated [8], i.e. either fully or effectively gaped. At the 
same time near the (tt, 0) momentum ARPES reveals a hump at ~ 100 meV 
[7, 5], which is consistent with large pseudogap scale T*. Spectral proper- 
ties for La2- x Sr x Cu204 (LSCO), as revealed by ARPES [6] appear to some 
extent different from BSCCO, presumably due to the crucial role of stripe 
structures in the LSCO in the regime of intermediate doping. Still they 
again reveal a truncated FS at low-doping and even the existence of QP 
along the nodal direction. 

The prototype single-band model relevant for cuprates which takes ex- 
plicitly into account strong correlations is the t-J model, derived originally 
by Chao, Spalek and Oles [9] 
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occupancy of sites, i.e., 

4 = (l-n,,_ s )4. (2) 

Longer range hopping appears to be important for the proper description of 
spectral function in cuprates, in particular it is invoked to explain the dif- 
ference between electron-doped and hole-doped materials both in the shape 
of the FS at optimum doping materials [10] as well as for the explanation 
of the ARPES of undoped insulators [11, 10], we consider besides Uj = t 
for the n.n. hopping also tij = if for the n.n.n. hopping on a square lattice. 
Note that t' < for the hole-doped cuprates. 

There have been so far numerous theoretical studies of the t-J model 
on square lattice and related Hubbard model at large Coulomb repulsion 
U S> t , as relevant to cuprates, using both analytical approaches as well 
as numerical techniques for finite size systems. Still analytical approxima- 
tions to spectral properties have proved to be very delicate, in particular 
with respect to the question of the emerging pseudogap at lower doping. So 
there are much fewer studies which give some answers on latter questions 
within microscopic models close to cuprates. The importance of AFM spin 
correlations for the emergence of the (large) pseudogap is found in the nu- 
merical studies [12, 13, 14] and in phenomenological model studies [15]. The 
renormalization group studies of the Hubbard model [16] also reveal the in- 
stability of the normal Fermi liquid close to the half-filled band (insulator) 
and a possible truncation of the Fermi surface. 

In the following we describe some evidence for the pseudogap within the 
t-J model obtained via finite-size studies and a novel approach to spectral 
functions using the method of equations of motion. 

2. Evidence for the pseudogap from numerical studies 

We introduced some years ago a numerical method, i.e. the finite- 
temperature Lanczos method (FTLM) [17, 14], which is particularly useful 
for studying finite-size model systems of correlated electrons at T > 0. The 
technical advantage of the method is that it is comparable in efficiency to 
ground state calculations. T > results for static and dynamical quantities 
are of interest in themselves, allowing to follow the T- variation of properties 
whereby some of them are meaningful only at T > 0, e.g. the entropy, the 
specific heat, the d.c. resistivity etc. On the other hand, the usage of finite 
but small T > represents a proper approach to more reliable ground state 
calculations in small systems. 

The most straightforward evidence for a pseudogap within the planar 
t-J model allowing the comparison with experiments appears in the uniform 
static spin-suceptibility Xo(T) [18, 14]. Results for various hole concentra- 
tions Ch = Nh/N in a system with N = 20 sites are presented in Fig. 2. 
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Fig. 2. Uniform susceptibility Xo(T) at several hole doping Ch- Results are for 
J/t = 0.3. 

It is evident that the maximum T* being related to the AFM exchange 
T* ~ 2 J/3 in an undoped AFM gradually shifts down with doping and 
finally disappears at 'critical' Ch = c* h ~ 0.15. Obtained results are qualita- 
tively as well as quantitatively consistent with experiments in cuprates, e.g. 



Another quantity relevant for comparison with the analytical approach 
further on is the single-particle density of states (DOS) M(lo) [13, 14]. We 
present in Fig. 3 numerical results for DOS [19] as a function of T for 
systems with N = 18, 20 sites and two lowest nonzero meaningful hole 
concentrations ~ 0.05,0.11. At smallest = 0.05 there is a pronounced 
pseudogap at uj ~ which closes with increasing T ~ T* ~ J. This 
again indicates the relation of this pseudogap with the AFM short-range 
correlations which dissolve for T > J. On the other hand, the pseudogap 
closes also on increasing doping since it becomes barely visible at Ch ~ 0.11. 

3. Spectral functions: Equation-of-motion approach 

In our analytical approach we analyze the electron Green's function 
directly for projected fermionic operators [20, 21], 



LSCO [3]. 
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Fig. 3. The DOS N{n + uj) at various T < J for hole doping: a) c h = 1/20 and b) 
c h = 2/18. 



which is equivalent to the usual propagator within the allowed basis states of 
the model. In the EQM method [22] one uses relations for general correlation 
functions 

u{{A;B)) u = ({A,B} + ) + {{[A,H];B)) U (4) 

applying them to the propagator G(k, lo) [20, 21] one can represent the 
latter in the form 

G(k,a;) = " — : r, (5) 

where a,Ck can be expressed in terms of commutators. It is important to 
notice that the renormalization a < 1 is already the consequence of the 
projected basis, 

« = ^E<fe>i}+> = ^( 1 + c ^ (6) 

while Ck represents the 'free' propagation emerging from EQM, 

Ck = ^{{[c ks ,H],cl s } + ) -( = -trnt-fr - 4.7^, (7) 
where rjj = a + (So-Sj)/a and 7k = (cos k x + cos k y )/2, 7^ = cos k x cos k y . 
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The central quantity for further consideration is the self energy 

E(k, W ) = ((C ks ;C+))L rr /«, iC ks = [c ks ,H] - ( k c ks , (8) 

and only the 'irreducible' part of the correlation function should be taken 
into account in the evaluation of S. EQM enter in the evaluation of ( k 
but even more important in C ks . We express the commutator in variables 
appropriate for a paramagnetic metallic state with (Sj) = 0, and we get 

[c ks ,H] = [(1 - |)e° - Jc e ]c ks + 

+ ~7Tf H m kq[s£qCk-q, s + 5^Ck_ q _ s - ^n q C k _q ;S ], (9) 
ViV q ^ 

where e k = —<itj k — 4t'j k is the bare band energy, hi = rii — c e and m k q is 
an effective spin-fermion coupling, 

™ k q = 2J 7q + e^_q. (10) 

One important achievement of the EQM approach is that it naturally leads 
to an effective coupling (not directly evident within the t-J model) between 
fermionic and spin degrees of freedom, which are essential for proper de- 
scription of low-energy physics in cuprates. Such a coupling is e.g. assumed 
in phenomenological models as the spin-fermion model [23, 15]. The essen- 
tial difference in our case is that m kc . is strongly dependent on k and q just 
in the vicinity of most relevant 'hot' spots. 



4. Self energy 

4-1. Undoped AFM and short range A FM fluctuations 

It is quite helpful observation that in the case of an undoped AFM 
our treatment of S and the spectral function reproduces quite successful 
SCBA equations [24] for the Green's function of a hole in an AFM. If we 
write EQM in the coordinate space for c\ s assuming the Neel state as the 
reference ni s = ±1, we get by considering only the t term, 

j n.n.i 

where we have also formally replaced dj s = Cj- s Sj. We note the similarity 
of Eq. (11) to the effective spin-holon coupling within the SCBA approach. 
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So we can follow the procedure of the evaluation of £AFM(k, to) within the 
SCBA in the linearized magnon theory [24] 

S AF M(k,w) = ^^M^ q G(k-q,u; + ^ q ), (12) 
q 

where w q is the magnon dispersion and Mk q = 4i(u q 7k- q + ^q7k) is the 
holon-magnon coupling which in general strong oc t, but vanishes near the 
AFM wavevector q = Q = (tt, tt). 

The advantage of the representation of EQM, Eq. (11), explicitly in 
spin and fermionic variables is that it allows the generalization to finite 
doping Ch > 0. We assume that spin fluctuations remain dominant at 
the AFM wavevector Q with the characteristic inverse AFM correlation 
length k = 1/£,afm- The latter seems to be the case for BSCCO as well as 
YB2Cu306+a;, but not for LSCO with pronounced stripe and spin-density 
structures [1]. For BSCCO and YBCO it is sensible to divide the spin fluc- 
tuations into two regimes with respect to q = q — Q: a) For q > k spin 
fluctuations are paramagnons, they are propagating like magnons and are 
transverse to the local AFM short-range spin, b) For q < k spin fluctu- 
ations are essentially not propagating modes but critically overdamped so 
deviations from the long range order are essential. 

At finite doping case we therefore generalize (at T = 0) Eq. (19) into 
the paramagnon contribution, 

£ pm (k,u,)=l]T [M^(k-q, W + Wq )+M^ qiq G + (k+q,a;- Wq )], (13) 

q,q>K 

where (^(k, uj) refer to the electron (uj > 0) and hole-part (uj < 0) of the 
propagator, respectively. We are dealing in Eq. (13) with a strong coupling 
theory due to t > w q and a selfconsistent calculation of S pm is required [21]. 
Also, resulting S pm (k, uj) and A(k, uj) as are at low doping quite asymmetric 
with respect to uj = 0, since G~ oc (1 — c/J/2 ~ 1/2 while G + oc c/j. 

4-2. Longitudinal spin fluctuations 

At Ch > the electronic system is in a paramagnetic state without an 
AFM long-range order and besides the paramagnon excitations also the 
coupling to longitudinal spin fluctuations become crucial. The latter re- 
store the spin rotation symmetry in a paramagnet and EQM (9) introduce 
such a spin-symmetric coupling. Within a simplest approximation that the 
dynamics of fermions and spins is independent, we get 

Eif(k,w) = -2^m k / / t/ wi,w 2 , 14 

a ^ H J J tt uj - uj 1 -uj 2 
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where g(co±, 0*2) = (l/2)[th(/3a;i/2)+cth(/3a;2/2)] and x is the dynamical spin 
susceptibility. Quite analogous treatment has been employed previously in 
the Hubbard model [25] and more recently within the spin-fermion model 
[15,26]. 

If we want to use the analogy with the spin-fermion Hamiltonian the 
effective coupling parameter mk q should satisfy mk, q = mk- q ,- q which is 
in general not the case with the form Eq. (10), therefore we use further on 
instead the symmetrized coupling 

m kq = 2J 7q + i(e k _ q + e k ). (15) 

In contrast to previous related studies of spin-fermion coupling [25, 15, 26], 
however, is strongly dependent on both q and k. It is essential that in 
the most sensitive parts of the FS, i.e., along the AFM zone boundary ('hot' 
line) where k = |Q — k|, the coupling is in fact quite modest determined 
solely by J and t' . Also, in the regime close to that of quasistatic x(q, w) 
the simplest and also quite satisfactory approximation is to insert for A the 
unrenormalized ^4° [15], the latter corresponding in our case to the spectral 
function without Sif but with £ = S pm . 

In the present theory spin susceptibility %(q, is taken as an input. 
The system is close to the AFM instability, so we assume spin fluctuations 
of the overdamped form [23] 

x"(q,w) oc — - — ^ — 2T 
(q 2 + k 2 )(uj 2 + uj 2 ) 

Nevertheless, the appearance of the pseudogap and the form of the FS are 
not strongly sensitive to the particular form of x"(q, u) at given character- 
istic k and lo k . 



(16) 



5. Pseudogap 

We first establish some characteristic features of the pseudogap and the 
development of the FS following a simplified analysis. We note that S pm 
induces a large incoherent component in the spectral functions at to <C 
and renormalizes the effective QP band relevant to the behavior at to ~ 
and at the FS [21]. It can also lead to a transition of a large FS into a small 
hole-pocket-like FS at very small Ch < c* h <^. 1. Nevertheless, the pseudogap 
can appear only via Sif. Therefore we here take into account S pm only via 
an effective band e k . The input spectral function for is thus 

4» f (k, w ) = <^ + /i - e ^). (17) 
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We restrict our discussion to T = and to the regime of intermediate (not 
too small) doping, where defines a large FS. The simplest case is the 
quasi-static and single-mode approximation (QSA) which is meaningful if 
w s «t,K«l, where we get 

The spectral function shows in this approximation two branches of E^, 
separated by the gap which opens along the AFM zone boundary k = \<lafm 
and the relevant (pseudo)gap scale is 

A£ G = |A kAFM | = ^-y/T s \2J - U'cos 2 k x \. (19) 

does not depend on t, but rather on smaller J and in particular t'. For 
t' < the gap is largest at (vr,0), consistent with experiments [7, 5]. 

QSA yileds a full gap corresponding, e.g., to the case of a long-range 
AFM state. Within the simplified effective band approach, Eq. (17), it is 
not difficult to evaluate numerically £if beyond the QSA, by taking ex- 
plicitly x"(q) w )> Eq. (16), with k > and uj k ~ 2Jn. For illustration, 
we present results characteristic for the development of spectral functions 
varying two most sensitive parameters k and fi, which both simulate the 
variation with doping, e.g. one can k take in accordance with experiments 
[1] and numerical results on the t-J model [27] as re ~ \fch- 

In Fig. 3 we present results for A(k, uj = 0) at T = for a broad range 
of k = 0.01 — 0.6. Curves in fact display the effective FS determined by 
the condition G _1 (kF,0) ~ 0. At the same time, intensities A{k, uj = 0) 
correspond to the QP weight Zp at the FS. At very small k = 0.01 we see 
a small (hole-pocket) FS which follows from the QSA, Eq. (18). Already 
k ~ 0.05 destroys the 'shadow' side of the pocket, i.e., the solution G^ 1 = 
on the latter side disappears. On the other hand, in the gap emerge now 
QP solutions with very weak Zp <C 1 which reconnect the FS into a large 
one. We are dealing nevertheless with effectively truncated FS with well 
developed arcs. The effect of larger k is essentially to increase Zp in the 
gapped region, in particular near (ir, 0). Finally, for large k = 0.6 corre- 
sponding in cuprates to the optimal doping or overdoping, Zp is essentially 
only weakly decreasing towards (ir, 0) and the FS is well pronounced and 
concave as naturally expected for t' < 0. 

We note in Fig. 4 that except at extreme k = 0.01 we get a large FS, 
whereby in the 'underdoped' regime k < k* ~ 0.5 the QP weight Zp is 
substantial only within pronounced 'arcs' and very small along the 'gapped' 
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Fig. 4. Contour plot of spectral functions A(k, u> = 0) at T = for various k in one 
quarter of the Brillouin zone. 



FS where Zp <C 1. Still, such a situation corresponds to a Fermi liquid, al- 
though a very strange one, where QP excitations exist everywhere along the 
FS and hence determine the low-energy properties of the 'normal' metallic 
state. 

It is quite remarkable to notice that in spite of Zp <C 1 the QP velocity 
vp is not diminished within the pseudogap. In fact it can be even enhanced, 
as seen in in Fig. 5 where the contour plot of vl(k, uS) is shown. Again, it is 
well evident in Fig. 5 that QP is well defined at the FS, while it becomes 
fuzzy at uj ^ merging with the solutions E^, respectively, away from the 
FS. The effect of large vf in the pseudogap, which is essential for the low-T 
thermodynamics, can be only explained with a crucial k dependence of S. 
Assuming that k enters only via e(k) we can express vf renormalization as 



K de'z^ 7 8e ^=o ~ w(uj k + w) ' l " 



While Zp/Z e f <C 1 in the same case, dT,'/de compensates or even leads to 
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an enhancement of vp. In the case u> K w <C A 2 we get 



r^i — 

„ef 



2J 



(21) 



Final vp is therefore not strongly renormalized, since 2 J and are of 
similar order. Furthermore, vp is enhanced in the parts of FS where is 
small, in particular near (ir, 0) point. The situation is thus very different 
from 'local' theories where S(k, to) ~ S(cj) and the QP velocity is governed 
only by Zp. In our case the 'nonlocal' character of S(k, oS) is essential in 
order to properly describe QP within the pseudogap region. 




(3^/8,7tt/8) 



Fig. 5. Contour plot of spectral functions A(k, ui = 0) across the FS in the pseudo- 
gap regime. 

^From A(k,u}) we can calculate the DOS M{uj) = (2/N) £k A(k, u). 
At the Fermi energy u ~ and low T one can express the DOS also as 
(assuming the existence of QP around the FS), 

The contribution will come mostly from FS arcs near the zone diagonal 
while the gapped regions near (vr,0) will contribute less due to Z(kp) <C 1. 
Results in Fig. 6(a) show the development of JV(0) with k. As expected the 
DOS decreases with decreasing k simulating the approach to an undoped 
AFM. The DOS is measured in cuprates via angle integrated photoemission 
spectroscopy, e.g. for LSCO in [28], as well as via the scanning tunneling 



Pseudogap printed on February 1, 2008 



13 



microscopy (STM) [29]. Our results are qualitatively consistent with these 
experiments (showing M(0) scaling with doping) as well as with numerical 
results on the t-J model [13, 14], as shown in Fig. 3. It is well possible 
that within photoemission experiments the matrix elements are essential, 
so we present in Fig. 6(a) also results weighted by the matrix element 
w(k) = (cosk x — cosky) 2 , which originates from the interplanar hopping 
as proposed for the c-axis conductivity [30]. Such weighted results show a 
stronger dependence on doping due to enhanced influence of the pseudogap 
region near (tt,0), which could be even closer to experimental findings [28]. 
In Fig. 6(b) we show also the average Z av along the FS, as well as the QP 
DOS, defined as 



The decrease Z av has similar dependence at the DOS M(0). It is however 
quite important to notice that smaller doping (decreasing re) leads also to a 
decrease of Mqp. This is consistent with the observation of the pseudogap 
also in the specific heat in cuprates [4], since Mqp oc 7 = Cy/T. We note 
here that such a behavior is not at all evident in the vicinity of a metal- 
insulator transition [1]. Namely, in a Fermi liquid with (nearly constant) 
large FS one can drive the metal-insulator transition by Z w — > 0. Within 
an assumption of a local character S(w) this would also lead to vp — ► and 
consequently via Eq. (23) to Mqp — > 00. Clearly, the essential difference in 
our case is that within the pseudogap regime S(k, u) is nonlocal, allowing 
for a large (not reduced or even enhanced) velocity within the pseudogap 
regime, hence a simultaneous decrease of M(0) and Mqp- 

It is also important to understand the role of finite T > 0. The most 
pronounced effect is on QP in the pseudogap part of the FS. The main 
conclusion is that weak QP peak with Zp <C 1 at T = 0, as seen e.g. in 
Fig. 5, is not just broadened but entirely disappears (becomes incoherent) 
already at very small T > T s <^ J, e.g. T s ~ 0.02 t for the situation in 
Fig. 5. This can explain the puzzle that ARPES experiments in fact do not 
observe any QP peak near (tt, 0) in the underdoped regime at T > T c [7, 5]. 



In this paper we have presented our results for spectral functions and the 
pseudogap within the t-J model, which is the prototype model for strongly 
correlated electrons and for superconducting cuprates in particular. Here we 
first comment on the validity of the model and results in a broader perspec- 
tive of strange metals and materials close to the metal-insulator transition. 
The physics of the t-J model at lower doping levels is determined by the 
interplay between the magnetic exchange (dominating the undoped AFM 




(23) 



6. Conclusions 
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Fig. 6. (a) DOS Af(0) and weighted DOS A/" w (0) vs. k. (b) Average QP weight Z, 
and QP DOS Af QP vs. k. 



insulator) and the itinerant kinetic energy of fermions (being dominant at 
least in the overdoped regime). Since itinerant fermions prefer a ferromag- 
netic state, the quantum state at the intermediate doping is frustrated, the 
quantum frustration showing up in large entropy, pronounced spin fluctu- 
ations, non-Fermi liquid effects etc. Evidently, this is one path towards 
the metal-insulator transition, but definitely not the only one possible. In 
this situation, fermionic and spin degrees of freedom are coupled but both 
active and relevant for low-energy properties. This is just the main con- 
tent and assumption of the presented theory for spectral functions and the 
pseudogap. 

The EQM approach to dynamical properties seem to be promising since 
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it can treat exactly the constraint which is essential for the physics of 
strongly correlated electrons. It has been recently also applied by present 
authors to the analysis of spin fluctuations and collective magnetic modes 
at low T. Our approximation for the self energy within the EQM approach 
deals with the normal paramagnetic state and treats the model as a coupled 
system (with derived effective coupling) of fermions with spin fluctuations, 
where close to the AFM ordered state both transverse and longitudinal spin 
fluctuations are important. Other contributions should be considered, e.g., 
the coupling to pairing fluctuations, in order to treat the superconducting 
state. 

In this paper we present only the results of the simplified pseudogap 
analysis. The results of full self-consistent treatment are qualitatively simi- 
lar [21]. Within the present theory the origin of the pseudogap feature is in 
the coupling to longitudinal spin fluctuations near the AFM wavevector Q 
which determine the QP properties in the 'hot' region, i.e. near the AFM 
zone boundary. The pseudogap opens predominantly in the same region 
and its extent is dependent on J and t' but not directly on t. Evidently the 
pseudogap bears a similarity to a d-wave-like dependence along the FS (for 
t' < 0) being largest near the (ir, 0) point. The strength of the pseudogap 
features depends mainly on k. It is important to note that apart from ex- 
tremely small k we are still dealing with a large FS. Still, at k < k* ~ 0.5 
parts of the FS near (ir/2, ir/2) remain well pronounced while the QP weight 
within the pseudogap part of the FS are strongly suppressed, in particular 
near zone corners (vr,0). 

The QP within the pseudogap have small weight Zp <C 1 but not dimin- 
ished (or even enhanced) vp, which is the effect of the nonlocal character of 
S(k, co). A consequence is that QP within the pseudogap contribute much 
less to QP DOS A/gp. This could be plausible explanation of a well known 
theoretical challenge that approaching the magnetic insulator both DOS, 
i.e. AA(0) and Aqp vanish. 

We presented results for T = 0, however the extention to T > is 
straightforward. Discussing only the effect on the pseudogap, we notice 
that it is mainly affected by k. So we can argue that the pseudogap should 
be observable for k(c^,T) < k* ~ 0.5. This effectively determines the 
crossover temperature T*(qJ. In the region of interest k is nearly linear in 
both T and Ch so we would get approximately 

T*~T *(1 -<%/<£), (24) 

where Tq ~ 0.6 J and c* h ~ 0.15. 

As described in previous sections, several features of our theory, regard- 
ing the development of spectral functions, Fermi surface and pseudogap, are 
at least qualitatively consistent with experimental results of normal state 
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properties in cuprates. However, further study within the present formalism 
is necessary in order to explore possible closer quantitative agreement with 
experiments as well as the emergence of superconductivity within the same 
model. 
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